Evaluating the partition function for systems with long range interactions 
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We express the partition function for an equilibrium system of interacting particles in the canonical 
ensemble as a functional integration over the particles' density field. We outline a method to evaluate 
the partition function by expanding around a saddle point. The saddle point is shown to be the 
solution of the equivalent mean-field theory. Leading corrections to the mean-field theory takes the 
form of a Gaussian integral. 
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The partition function for a system of interacting par- 
ticles in equiUbriuni can be written down quite easily. 
Using it to evaluate physically relevant quantities, how- 
ever, generally requires aggressive approximations. For 
example, expansion in powers of the density is only good 
for short-range interactions and low density El . Few ap- 
proximation schemes are appropriate for strongly inter- 
acting systems characterized by long-ranged forces, such 
as Coulomb interactions and gravitational forces. 

The method most often used to address systems with 
long range interactions is the mean-field theory, in which, 
for any chosen particle, the effects of all other particles 
are replaced by a mean-field which is not affected by this 
particle. Its advantage lies in its simplicity and physical 
clarity. Due to the long-range nature of the interactions, 
mean-field theory is often a good approximation. How- 
ever, its inadequacy is shown in experiments ^-^■ 

The lack of thermodynamical self-consistency of the 
mean-field theory was pointed out nearly seventy years 
ago [HJ^ . Efforts were then made to develop rigorous ap- 
proaches for treating such systems by Onsager la]. Kirk- 
wood Ig], and Bogolyubov ||l[. These methods, from dif- 
ferent mathematical view points, however, all relied on 
detailed descriptions of the microscopic states. Attempts 
to treat such systems by finding descriptions involving 
increasing number of particles are undermined by the 
notorious closure problem. Consequently these methods 
failed to make any practical predictions, nor an estima- 
tion for the accuracy of the mean-field theory. 

Philosophically, the failure of these methods is due to 
the fact that a system with long range interactions is 
highly coupled — any part of it is affected by the changes 
in the whole system — and cannot be described effectively 
by schemes focusing on local effects. This suggests an 
approach which treats such a system as a whole. The 
relevancy of the system as a whole also determines that 
the detailed local information becomes more and more 
irrelevant with increasing system sizes, and thus can be 
replaced for a big system through coarse-graining. 

This Letter proposes a new approach based on these 
considerations. We express the partition function as a 
functional integration over the particles' density field. 
The saddle point of this expression turns out to be equiv- 



alent to the mean-field theory. Corrections to the mean- 
field theory then can be expressed as an expansion around 
this saddle point. 

Let us consider a system with N classical particles, 
volume V and temperature T. Its partition function in 
the canonical ensemble is. 
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where U{qi ■ ■ ■ q^) is the potential energy due to external 
fields or interactions among particles. In the following, 
we will concentrate on the integral part in the above ex- 
pression, so let us write. 



Zs = 



-'^^dqi---dqN. 
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Equation (y) is written in terms of integrations over 
the particles' positions. We want to instead express Zg 
in terms of functional integrations over the particles' den- 
sity field. To derive such an expression, let us start with 
the ideal gas — the factor exp(— /?[/) can be simply in- 
serted into the result later. 

For an ideal gas, Zg = V^ . Let us divide V into 
M cells, with volumes Vi,V2,- ■ ■ ,Vm- A state n is 
represented by the number of particles in these cells, 
71 = (rti, 712, • • • , «A/)- The division must satisfy the con- 
dition that there are many particles in each cell, i.e. rii 
is big. The probability weight of such a state is denoted 
by W{n) and it can be obtained through the following 
considerations. 
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where the / above the summation sign indicates that par- 
ticle number is preserved: '^rii ~ N . 



We can evaluate Eq. Q using Stirling's approxima- 
tion, 



because rii is big. Then 
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i=l 

And we express this in terms of the density pi = rii/Vi 
of particles in cell i, 
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We use this form of W{n) to evaluate Z^, and change 
from summations over different configurations of rii to 
integrations over continuous variables ni, 
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For general cases, we simply insert the factor e '^ , so 
Zs = \/2ttNN^ I Dpe-'^f, (14) 



where 



Dp = 
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f = U + kT plogpdV. 
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Equations (|l4|-|l6|) are the desired expression for Zg and 
are the central result of this letter. 

When we replace the position variables of the parti- 
cles with the coarse-grained density distribution, certain 
information is lost. In particular, U is calculated not 
for the actual positions of the particles, but instead for 
the coarse-grained density distribution. For this to be 
a good approximation, certain quantities such as field 
potential must not vary appreciably within a cell; each 
cell is treated as a uniform subsystem. This is an addi- 
tional condition to the one mentioned above that there 



are many particles in each cell. Finding an appropriate 
division of a macroscopic system which satisfies both con- 
ditions only is possible if the range of the interaction is 
long enough. 

To calculate Z^ by using Eq. (n4h directly is as difficult 
as by using Eq. (Q). However, the form of Eq. (|lj) sug- 
gests a different approach — first, find the saddle point, 
i.e. the minimum of /, and then expand around this sad- 
dle point. 

Let us assume po is the density distribution which min- 
imizes /. The corresponding minimum is /q. Then, 
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To study the saddle point, let us consider an example 
of an system of particles with Coulomb interaction, en- 
compassed by a surface E with surface charge a. For sim- 
plicity, let us assume all particles carry charge e. Then 
the interaction potential energy for a configuration p of 
particles is 



U^l f pe(t)dV + \ I (j(j)dS, 
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where (j) is the electric field potential, and is connected 
to p by the Poisson equation. 
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where n is the unit vector normal to the surface, and 
pointing outside the system. So, keeping in mind that 
/ 5pdV = 0, 
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The second term in the last expression vanishes because 
there is no charges outside of the system. Then the van- 
ishing of the first term leads to the Boltzmann distribu- 
tion. 



Po = Qoe 



(24) 



where the constant Qq ensures J podV — N. Now, more- 
over, we find that the saddle point corresponds to the 
Poisson-Boltzmann (PB) equation. 
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Even though the above example focused on the special 
case of Coulomb interactions, Eq. (|lq ) is completely gen- 
eral so that the saddle point always correspond to the 
system's mean-field theory. This is reasonable because 
physically, when we consider only the optimal density 
distribution and ignore all fiuctuations, we get the mean- 
field theory. Of course, the mean-field theory, though of- 
ten a good approximation, is neither complete, nor ther- 
modynamically self-consistent. In particular it ignores 
fiuctuations whose effects we will consider next. 

We account for fluctuations by expanding around the 
saddle point, p = po + 5 p. The first-order change in / 
vanishes, let us write the second order change in / as 
5^ f , and in U as 5^U . 6^f is a quadratic functional of 
6p, 
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and thus is always positive. For Coulomb and gravita- 
tional interactions, there are no higher order corrections 
to U than S^U. 

Now we have, from Eq. ([l4|). 
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In obtaining Eq. (p7|), we neglect terms to the order of 
{5p)^ in comparison to terms to the order of {5pY in the 
expansion of p log p, also in (Eq) , we write poi instead of 
Pi in the denominator. Both approximations are justified 
by the fact that each cell contains many particles. The 
prime above the integration sign in Eq. (|2^) indicates 
that particle number is preserved, J pdV — N . 
Now we can write. 
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Thus, we separate Zs into two parts: Zq describes the 
mean-field theory, while Z^ accounts for fluctuations. 



For non- interacting particles, 5^U = and the Gaussian 
integral for Zm gives Zm = 1. The mean field theory is 
exact for this case, as it should be. However, for inter- 
acting particles, Zm yields a nontrivial modification. 

For example, ref. investigates a model colloidal sys- 
tem of a charged sphere surrounded by counter-ions in- 
side a spherical shell, mimicking the Wigner-Seitz cell 
in a colloidal crystal. The mean-field (PB) theory pre- 
dicts that the free energy is a monotonicly decreasing 
function of the radius of the shell. However, Z^ con- 
tributes a increasing component to the free energy thus 
leads to a local minimum in the free energy's depen- 
dence on the shell's radius. This would imply an effective 
attraction between like-charge colloidal particles, which 
was observed in experiments P,Hl. 

If, instead of charged ions, we have neutral particles 
which interact with attractive gravitational forces in the 
above model, the mathematical expressions would be 
nearly identical except that the interaction is with a dif- 
ferent sign. Then immediately we see Z^ would lead to 
an extra effective repulsion. This may be relevant to the 
observation of anomalous gravitational effects |g] . 

Let us consider the effects of dividing the system into 
cells from a different view point. In order to derive Equa- 
tions (|29|-pl|) , we need the condition that in each cell there 
are many particles. Of course, the final results should not 
depend on the way how we divide the system. Essentially, 
we are considering the effects of density fluctuations. If 
we decompose these fiuctuations into modes with differ- 
ent wavelengths, then the size of the cell corresponds to 
a short-wavelength cutoff. When we calculate Eq. (|3l]), 
there are two contributions: one from d^U and the other 
from the integration of ^^^^ . The latter is totally insensi- 
tive to the way the cells are made — if we can ignore 6^U, 
we always get Zm = 1 no matter how we divide the sys- 
tem. For long range interactions, we would indeed expect 
S'^U to be very insensitive to short wavelength fluctua- 
tions. So the final result should be independent of the 
way the system is divided. 

Equations ([29|-pl|) remind us of a basic concept: The 
density distribution is not one of the macroscopic quan- 
tities such as temperature and volume which determine 
the equilibrium state of the system. It cannot be preset, 
but rather results from an ensemble average. Thus when 
there are interactions between particles, thermodynamic 
properties cannot be calculated from a fixed optimal den- 
sity distribution in a self-consistent way. The formalism 
developed here remedies this deficiency. 

The method developed here bears a passing resem- 
blance to the density functional theory (DFT) p-fl^] in 
that both approximations use density as the basic vari- 
able. However, DFT calculations typically are based on 
an optimal density distribution, since the goal is often to 
calculate a system's ground state. Here, density fluctua- 
tions are a essential component as they lead to Eq. (pl|). 
For finite temperature, DFT calculations based on one 



optimal density distribution cannot be thermodynami- 
cally self-consistent for reasons stated above. In that 
case, the method developed here may be helpful to ex- 
tend DFT calculations. 

Also, Ornstein derived the van der Waals equation in 
a very similar way |13|| . But again, the purpose there was 
to find the optimal density distribution and no fluctua- 
tions were considered. 

Netz and Orland recently formulated a field theoreti- 
cal approach ||Ij| to the same class of problems addressed 
herein. They obtained the PB equation as a saddle point 
and calculated corrections by an diagrammatic expan- 
sion. This approach is well developed in field theory. 
However, the convergence of the expansion relies on a 
small parameter which is usually missing when the in- 
teraction is long ranged. Actually, Coalson and Duncan 
had shown that the subsequent terms in such an expan- 
sion are of the same order |g5[ . The method developed 
here, on the other hand, exploits the long range of the 
interactions to coarse grain the system effectively. 

This work is motivated by the experimental results 
showing the inadequacy of the mean-field theory [^,||. 
For an application of the method described here to 
a model colloidal system, see 1^. The author wants 
to thank David Grier, Susan Coppersmith and Daniel 
Chung for helpful discussions. This work was supported 
in part by the MRSEC Program of the NSF under award 
number #DMR-9400379. 
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